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SOLUTION OF A SECOND-ORDER LINEAR SYSTEM 


BY MATCHED ASYMPTOTIC EXPANSIONS 
Mark D . Ardema 
NASA Ames Research Center 

SUMMARY 

For the purpose of giving a clear exposition of the method, matched asymptotic 
expansions (MAE) are used to obtain a first-order approximation to the solution of 
a singularly perturbed second-order system. A special case is considered in which 
the uniform asymptotic solution obtained by MAE is shown to converge to the exact 
solution. Ways in which the method can be used to solve higher— order linear systems, 
including those which are not singularly perturbed, are also discussed. 

INTRODUCTION 


One of the simplest singular-perturbation problems is that of finding an 
asymptotic solution to a constant-coefficient linear system, with one slow and one 
fast variable subject to initial conditions. Such a system is given by 


= Ax + By ; x(e,0) = a 
e ^ = Cx + Dy ; y(e,0) = 3 


( 1 ) 


where x and y are called the slow and fast variables, respectively, and e is 
a "small" parameter. In this paper we will obtain an asymptotic solution to this 
problem to first order in e on 0 < t ^ T < 00 by the method of matched asymptotic 
expansions (MAE) under the following two assumptions: 


( 2 ) 

(3) 

The first of these assumptions assures the existence of a solution and the second 
assures boundary-layer stability, the significance of which will become apparent 
later. 


E = AD - BC 0 
D < 0 



The purpose of solving such a simple example is to make as transparent as 
possible the use of the MAE method. This method was developed in connection with 
certain fluid mechanics problems and is discussed in depth in references 1 and 2. 
Reference 3 adapts and applies MAE to a wide range of problems involving systems of 
ordinary differential equations. 

There are four steps in the MAE method: (1) formulation of the outer problem 

and its asymptotic solution to the desired order; (2) formulation of the inner 
(boundary-layer) problem and its asymptotic solution to the desired order; (3) match- 
ing of the outer and inner solutions to obtain unknown constants of integration and 
the common parts; and (4) formulation of a composite solution that is uniformly 
valid in all dependent variables for 0 < t < T < °°. 

Although the MAE method is somewhat more cumbersome than other singular- 
perturbation techniques that have been developed in recent years, it has two 
advantages over the more recent methods. The first is that it explicitly splits the 
problem into two (or more) separate ones. The solutions of these problems may be 
useful in their own right, in addition to being the constituents of composite 
solutions. Further, the inner problem(s) satisfies the boundary conditions as 
opposed to most other methods. The second advantage lies in the great flexibility 
that is afforded in the form of the assumed stretching transformations and inner and 
outer expansions, as well as in the manner in which these expansions are combined to 
give the uniform composite solution. This flexibility is often desirable and 
occasionally necessary. 

The solution to equation (1) will be obtained to first-order, since determination 
of the zero-order solution is essentially trivial. The zero-order uniform solution 
consists of the solution to the reduced problem [Eq. (1), with e = 0 and the initial 
condition on y discarded] added to the solution of the zero-order inner problem 
[eq. (1), with t transformed to t - t/e and e = 0] . Steps three and four of the 
MAE method are then not needed. Further, the first-order solution gives much better 
numerical results than the zero-order one, unless e is ,f very small. 11 This point 
is dealt with in the final section ( Discussion and Generalizations ) . 
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By definition, f(e) is an nth-order asymptotic approximation to F(g) if 

iim | F(e) - f(e)| < 

e-K) n+1 

e 


OUTER SOLUTION 

The outer system associated with equation (1) describes the solution away from 
the initial point. It is simply equation (1) without the initial conditions: 


dx 
dt 

iz = 

dt 


Ax° + By° 


(4) 


Cx° + Dy° 


To solve this to first-order we assume that power-series expansions are valid and 
set 


O. . O , . , o, . 

x (e,t) = x (t) + x,(t)s 
o i 


o , N o , . . o , . 

y (e,t) = y (t) + y (t)c 
o 1 


in equations (4) and retain only first-order terms to get 


• ° j o 

dx dx_ 

1 « ° * » o _ o - _ o 

~rr" + ~T~ e = Ax + Ax_ e + By + By- c 
dt dt o 1 o 1 


dv o o o o 

e = Cx + Cx n g + Dy + Dy. g 
dt o 1 -^o J 1 


Equating the coefficients of e gives the zero-order outer problem 


(5) 


( 6 ) 


dx 

° A o , _ o 
—r— = Ax + By 
dt o J o 


0 = Cx° + Dy° 
o o 


(7) 
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The solutions to equations (7) under assumptions (2) and (3), are 


It 

x° = C° e D 
o o 


1 1 

o -0 C D 

y = -C — e 
o o D 


( 8 ) 


where C q is a constant of integration that is as yet undetermined. Since 
equations (7) are a first-order system, both of the initial conditions cannot be 
satisfied (unless, of course, 6 - -a — ) . Thus, at best, equations (8) are zero-order 
approximations only on 0 < t < T < °°, and no initial conditions have been imposed 
on the outer problem. 

Similarly, the first-order problem is obtained from the coefficients of e in 
equations (6) as 


^i 

dt 


^o 

dt 


. o , „ o 
Ax l + By l 


o o 

Cx 1 + Dy x 


(9) 


Since y°(t) is a known function from the solution of the previous order, this is a 
o 

first-order problem. Using (8), 


which has solution 


dx - 

1 E o _o CBE D 

— — = - x, - C — — e 
dt D 1 o 3 


o „o CE D C Co 

H ' - c o ^ e - 5 x i 




,0 CBE 


t e 


h 


( 10 ) 
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( 11 ) 


y l = 


■C° e“ * - C° § e” C + C° ^ t j> ' 
0 D 3 1 D ° D 4 


where C° 


is another as yet unknown constant of integration. 


INNER SOLUTION 


Because the purpose of the inner system is to model the solution near the initial 
point, its solution is required to satisfy both initial conditions. Postulating that 
the inner motion occurs on a time scale of order e, we stretch the independent 
variable by 


_t 

£ 


( 12 ) 


Substituting equation (12) in (1) gives the inner (boundary-layer) equations 


dx i . _ i. i, 

= e(Ax + By ) ; x (e,0) = a 


- Cx 1 + Dv 1 

QT 


Dy 1 ; y i (e,0) = 3 


(13) 


To solve this to first-order, assume a power-series expansion and set 


x i (e,x) = x i (x) + x^(x)e 
o 1 


y i (e,t) = y*(x) + y*(x)e 


in (13) to get 


(14) 


dx 1 

dx* 

-{ -i 

U . X 

dx dx G 

= Ax e + By e ; 

o J o 

i 

i 


dy o 

d y, 

1.1 i i 

— + -T*- e 
dx dx 

" CXjj + Cx 1 e + Dy Q + Dy 1 e ; 


; x*(0) + x*(0)e = 

o 1 


= 8 


(15) 
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The zero-order problem from (15) is 


0 ; x i (0) = ct 

o 

K + d 4 i ^<°> - s 

with the solution 


i 

x = a 
o 



dx 
o 

dt 


^o 

dx 


( 16 ) 


(17) 


Note that x 1 is a constant. In effect, the boundary-layer motion is so rapid that 
to zero-order the slow variable has not had time to begin its motion. 


The first-order problem is 


^i 

dx 


Ax 1 + By 1 
o o 


x*(o) 


= o 


dy, 

dx 


CxJ + Dy^ ; y^(0) 


= 0 


( 18 ) 


which has the solution 



Note that to first-order there is a variation in the slow variable x in the 
boundary layer. 
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MATCHING 


Matching is the key step in the MAE method. It identifies any faulty assumptions 
regarding the form of the asymptotic expansions and stretching transformations and 
serves to determine the unknown constants C° and C°. Matching essentially requires 
that the behavior of the outer solution as t -*■ 0 is the same as that of the inner 
solution as t ®, that is, the outer solution extended into the inner region must 
agree with the inner solution extended into the outer region. This "limit matching 
principle" may be stated as 

lim [x°(e,t) - x^(e,t)] = 0 (20a) 

t-*0 

x->°° 

e -*0 

lim [y°(£,t) - y 1 (£,T)] = 0 (20b) 

t-vO 

X~HX> 

e-K) 


These equations are to be regarded as shorthand for the requirement that the outer 
and inner solutions must agree in an "overlap region" between the inner and outer 
regions. 


The limiting behaviors of x° and y° are obtained by expanding equations (8) 
and (10) about t = 0; for small t, 


o 

o 

c i 

o 


+ D c + • • •) 

C o D (- 1 + D C + ' ' ') 

+ (-•••) 


t + 


( 21 ) 


• •) + c d “ (, * ...j 


For large t, (17) and (19) are approximately 
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i 

x % a 
o 


0 

1 


C 

' D “ 


I ax _ ! /g + CcA 
D D y D J 


( 22 ) 


i ECa ECa . 

y i * ‘ TF T - 3" + 


BC ( Ca\ 

D 2 V D ) 


Dx 

By assumption (3), terms involving e do not appear in (22) 


Using equations (5), (14), (21), and (22) in equation (20a), we now match x to 
first-order: 


im |c° + C° ^ 
+0 ( ° ° D 


lim 

t->0 

E +0 


— t + ... + e 


ft + f 


t +' 


C°^t- 
O _ J 


D~ 


- a - e 


E 

-CT 


') 

' 5 (‘ * t)]J ' " 


(23) 


Noting that c goes to zero faster than t, and using equation (12), equation (23) 
can be true only if 



o 

O 0 

i 

p 

n 

o 

(coefficients 

of 

O Ov 

t e ) 

(24a) 

o 
o o 

E 

D 

E 

- a = 0 

(coefficients 

of 

.1 o 
t e ) 

(24b) 

*5 

( B + 

O 

II 

IP 

(coefficients 

of 

. ° lv 
t e ) 

(24c) 


Equations (24a) and (24c) give 


C° = a 
o 




(25) 
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and (24b) is then satisfied trivially. Terms in equation (23) that involve higher 
powers of t m e n than t°e 1 are not matched at this order of approximation. 

It is now clear why assumption (3) is necessary. If D > 0, then y^ grows 
exponentially and cannot be matched by any algebraic term. It is in fact a general 
requirement of singular-perturbation analysis that the zero-order boundary-layer 
equations, equations (16) in our case, be asymptotically stable. In the MAE method, 
this requirement arises naturally from the matching relation. 

The matching relation for y to first-order is 



This implies 


o CE o C ECa 

0 D 3 1 D D 3 



0 

0 

0 


(27) 


which agrees with (25). In general, matching the slow variables is sufficient to 
determine all unknown constants of integration. 


The fact that we have been able to match both x and y validates our 
selection of the linear stretching transformation (12) and power-series expansions (5) 
and (14). If matching had not been possible with these assumptions, more general 
transformations and expansions would be required. 
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From equations (25), (8), and (10) the outer solution is 



FORMATION OF COMPOSITE SOLUTION 


( 28 ) 


We now have a representation for the solution of (1) near t = 0, as given by 
equations (14), with (17) and (19), and a representation away from t = 0, as given 
by (5), with (28). However, it is generally more useful to have one solution that 
is valid everywhere, and this is the purpose of composite solutions. The most common 
composite solution (in principle the number of solutions is infinite) is -the additive 
one formed by simply adding the inner and outer solutions. The result must be 
adjusted by subtracting out the "common part," that is, the portion of the solutions 
that explicitly cancels in the matching relation; otherwise, this portion would be 
added in twice. In particular, the initial conditions would not be met. Thus, the 
additive composite has the form 


To first-order. 


x a (e, t) = x°(c,t) + x 


y a (e,t) = y°(e,t) + y 


l H)- 


CP x (e,t) 


CP (e,t) 


(29) 


a o . o , i , i __ 
x, = x + x.e + x + x. e - CP 

1 O 1 O 1 X, 


a o . o , i . i __ 
y l = y o + y l e + y o + y l E ' CP y 


(30) 
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From equations (23) and (26), or equivalently (24) and (27), the common parts 
to first-order for the problem at hand are 


cp x 1 = “ + !“ t -!( 6 + f) E 

„„ C EC . EC BC ( 0 . Ca \ 

c V-s 0 -7 at -7“V( #t ' 5 7 


(31) 


Note that although the matching of y is not required for determining constants, 
it is required for determining common parts. According to (30), the first-order 
additive composite solution is 


x, = 


a e - e 


+ 


E D 

U /■ ^ C«\ ± CBE I D * B / Ca \ 1 

I D + “dV + “ “^3“ t J e +e S( 6 + Tj e 

— r 2 *i - 1 

a £ e »\ £ [_^ + M( 6 + ^ + a ^ t p 


(32) 


(33) 


This solution is a first-order, uniform asymptotic approximation to the problem (1) , 
that is, if x(e,t) and y(e,t) are the exact solutions, then 


x(e,t) - x (c,t) 

lim 5 < 00 

e-*-0 e 


y(e,t) - y (e,t) 

lim 2 < 00 

e+0 e 


on the interval 0 ^ t < T < 00 . 
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A SPECIAL CASE 


In this section, a special case is considered for the purpose of easily 
determining higher-order approximations and comparing the approximate solution with 
the exact one. Specifically, consider (1) with C = 0: 


^ ^ + Sy 


x(e,0) = a 


iz = 

dt 


Dy 


y(£,0) = 6 


(34) 


In this system, the second equation is uncoupled from the first and the exact solution 
is easily obtained as 


fit' 

e At BB | e C - e G 

x = a e - 


(D/e - A) 


e e 

y = 3 e 


Solving the outer and inner problems associated with (34) gives 


(33) 


o 

X. 

J 


C° e At 


j = 0,1,2, 


i . , BB / Dt A 

X 1 = A«t + — ^e - 1 J 

i 1 ,2 2 ABB . ABB / Dt A 

x 2 = 2 A ° T - “5“ T + -J \ e ~ X ) 


(36) 


y = B e 


Dt 


i 

y j 


= 0 


» j 1 > 2 » • 
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Forming the additive composite for y gives 


so that 


St 


j - 0,1,2, 


a e 

y i = y 


5 j 0,1,2, • • 


(37) 


(38) 


This will happen whenever the fast equation does not contain slow variables and the 
initial condition on y does not depend on e. 


Forming the additive composite for x to second-order. 



From this we conjecture that 



This clearly shows the need for f e (A/D) | < 1. By the binomial expansion: 


(39) 


(AO) 


1 

(D/e - A) 



(41) 


From (35) and (41), 



( 42 ) 
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Comparing equations (40) and (42) we see that 


lim xf = x 6 (43) 

so that the uniform asymptotic expansion resulting from the method of MAE converges 
to the exact solution. 


DISCUSSION AND GENERALIZATIONS 

We now return to the problem stated in the Introduction. The fact that the 
solution (32) obtained by MAE is a first-order, uniform, asymptotic approximation 
to the exact solution gives no assurance either that the asymptotic expansion is 
convergent or that it is numerically close to the exact solution. The first question, 
that of convergence, is relatively unimportant. (This is somewhat surprising in view 
of the heavy emphasis placed on convergence in elementary mathematics.) The second 
question, that of numerical accuracy, is obviously of great importance in practical 
applications and will be taken up now. 

The solution (32) consists of two types of terms, those with factor exp[(E/D)t] 
and those with factor exp[(D/e)t]. The first type arises from the outer solution 
and the second from the inner. The more rapidly the inner terms decay relative to 
the outer terms the better the approximation (32) will be. Thus, it is not the 
absolute magnitude of e that is important for numerical accuracy but rather the 
size of | D/e | relative to | E/D | . The larger |D/e| is relative to |e/d|, the 
better will be the approximation. We write this accuracy requirement as 


E 

D 



(44) 


Now suppose e = 1, that is, the system is given by 


dx 

dt 

Ax + By ; 

x(e,0) = a 

(45a) 

ix = 

dt 

Cx + Dy ; 

y(e,0) = B 

(45b) 
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Then (44) becomes, using (2), 


AD - BC 


<D 4 


(46) 


This will be true if 


|D| > Max(|A| , | B | , | C | ) (47) 

Thus we can solve the nonsingularly perturbed problem (45) by MAE, provided (47) 
holds. The solution to first-order is (32) with e = 1. The condition (47) may 
also be deduced by studying the eigenvalues of the coefficient matrix of (45). 

These eigenvalues are 


V X 2 = 1 { A + D ± [ (A " D)2 + 4BC] J (48) 

If (47) is satisfied, one of these eigenvalues will be approximately equal to D, 
a large negative number, and the other will be relatively small in magnitude, 
exactly the situation we require. 


Solution of (45) by MAE may be viewed another way. Let 

Max ( | A | , 1 B | , [Cl) 

l D l 

Multiply (45b) by e to get the system 

“ Ax + By ; x(e , 0) = a 

e ^ = C’x + D'y ; y(e,0) = $ 


(49) 


(50) 


and let 


E' = AD' - BC' 


(51) 



where 


C» = e c ; D’ = eD (52) 

We now have a problem to which the MAE method may be applied, and the first-order 
solution is given by (32), with C, D, and E replaced by C’, D', and E' . 
Transforming back into unprimed quantities then gives (32), with e = 1. Therefore, 
the same answer is obtained by simply artificially inserting e in front of dy/dt 
in (45b), applying MAE, and then setting e = 1. This technique greatly broadens 
the applicability of singular-perturbation methods. 

Since there is no distinction between x and y in equations (45) , we conclude 
that for a valid and useful solution of (45) by MAE (or any other singular-perturbation 
method) it is sufficient that either 


D | ^ Max (| a|,|b|,|c|), and D < 0 


or 


A | ^ Max ( | D | , | B | , |c|), and A < 0 


(53) 


(54) 


It is obvious that these conclusions may be easily extended to higher-order 
systems. Consider 


dx 

— = A(e,t)x + B(e,t)y ; x(e,0) = a 

dt = 

dy 

— = C(e,t)x + D(e,t)y ; y(e,0) = B 

dt 


(55) 


oo 

where A(e,t), B(e,t), C(e,t), and D(e,t) are of class C in t and have 
asymptotic power-series expansions in e. A valid and useful solution of this 
problem can be obtained by MAE, provided the solution consists of combinations of 
fast decaying and slow modes. This will occur when the eigenvalues of D 
(or equivalently, of A) have real parts that are relatively large (in absolute value) 
negative numbers. Therefore, we can give the following general procedure for solving 
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linear systems by MAE: (1) put the system in a form such that eigenvalues of D 
have the required properties by rearranging the equations or transforming the 
variables, or both; (2) insert e in front of djr/dt to create a singularly 
perturbed system; (3) apply the four steps of MAE; and (4) set e = 1. 
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